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ABSTRACT 


A one  dimensional  time-marching  model  to  solve  the  dynamic 
behavior  of  a compression  system  is  presented.  The  computational 
meLhod  is  described  in  some  details  and  it  is  applied  to  a three-stage 
axial  compressor  with  different  plenum  configurations.  This 
computational  model  for  the  simulation  of  a multistage  axial 
compressor  is  used  to  investigate  the  insertion  of  an  active  control 
system. 
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duct  section,  sound  speed 
internal  energy 

mass,  momentum  and  energy  flux  vector 
blade  force,  external  force  due  to  friction 
enthalpy,  enthalpy  of  indeed 
specific  heat  ratio  Cp/cv 
duct  length,  compressor  length 
number  of  computational  sections 
mass  flow  bleed 
static  pressure 

conservative  flow  variables  vector 
gas  constant 
time,  static  temperature 
velocity,  peripheral  speed 
Cartesian  co-ordinates 
plenum  volume 
work  exchange 

forward  and  backward  finite  differences 
flow  coefficient,  external  heat  flux 
pressure  coefficient 
valve  coefficient 
density 

total  quantity 
inlet,  outlet  sections 


INTRODUCTION 

The  need  for  a wider  operability  range  of  modem  compressors  for 
gas  turbine  applications  is  the  "prime  mover"  for  a large  research 
effort  that  is  nowadays  undertaken  in  many  international  laboratories, 
organisations,  companies  and  universities.  This  effort  is  both 
experimental  and  computational.  The  overall  behaviour  of  a gas 
turbine  compressor  outside  its  stability  range  has  been  understood  and 
can  be  modelled  using  a simplified  lumped  parameters  approach 
(Greitzer,  1976,  Baghdadi  et  al.,  1982).  Many  detailed  experimental 
analysis  are  nowadays  focused  on  stall  inception  (Day  et  al,  1999, 
Camp-Day  1998,  Spakovsky  et  al.,  1999)  to  understand  the  fluid- 
dynamic  mechanism  of  formation  and  to  be  able  to  improve  its 
computational  modelling.  On  the  other  hand  the  computational  models 
for  the  analysis  of  a multistage  compressor  in  steady  flow  (through- 
flow  or  3D  Navier- Stokes)  are  not  directly  applicable  to  unsteady 
(transient  or  dynamic)  flows  because  they  are  either  inappropriate 
(standard  through-flow)  or  they  require  an  excessive  computational 
effort  (3D  Navier-Stokes).  The  most  simplified  approach  is  the  zero 
dimensional  lumped  parameter  technique  that  considers  each 
component  of  the  compression  system  as  a node  and  by  writing  the 
balance  for  mass  flowr,  momentum  and  energy  it  results  in  a set  of 
differential  equations  to  be  solved  with  respect  to  time  (Greitzer, 
1976,  Massardo  et  al.  1989,  Botros,  1994). 

The  one  dimensional  model  for  the  analysis  of  the  unsteady  flow  in  a 
compression  system  can  be  a good  compromise  between  the  accuracy 
in  capturing  the  main  system  performances  and  the  computational 
effort.  This  approach  introduces  the  conservation  equations 
(continuity,  momenLum  and  energy)  for  a continuum  and,  after  a 
domain  discretization,  are  integrated  using  a time-marching  technique; 
the  effect  of  blades,  mass  bleeds,  friction  etc....  are  introduced  as 
external  body  forces.  This  technique  has  been  considered  and  a time- 
marching technique,  previously  developed  for  2D/3D  turbomachinery 
flows  (Cravero,  1995),  has  been  converted  in  ID  form  with  the 
insertion  of  the  appropriate  external  forces  to  model  the  dynamic  of 
the  compression  system.  The  ID  lime-marching  approach  is  preferred 
over  the  lumped  parameter  technique,  because  it  allows  the  analysis 
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inside  the  compressor  (interstage  analysis)  and  this  is  an  important 
feature  when  an  active  control  strategy  must  be  defined.  In  fact  the 
control  devices  and  the  transducers  should  be  placed  in  the  most 
efficient  way  (Escuret-Elder,  1993,  Gysling-Greitzer,  1994,  Montazeri 
et  al.,  1996,  Sun-Elder,  1998)  and  this  can  be  done  if  a reliable 
technique  for  the  compression  system  analysis  is  implemented  Active 
control  of  compressor  instabilities  (surge)  has  received  attention 
recently  (Epstein  et  al.,  1989);  it  aims  to  allow  the  compressor  to 
operate  beyond  the  surge  line  suppressing  the  flow  instabilities  by 
means  of  a feedback  control  system.  The  design  of  active  control 
systems  for  multistage  axial  compressors  is  still  a challenging  topic;  it 
requires  a good  representation  of  the  compression  system  under 
unsteady  flows  in  order  to  detect  the  instability  and  to  help  the 
optimisation  of  the  control  system  (transducer  and  actuator 
positioning). 

COMPUTATIONAL  MODEL 


The  geometrical  domain  is  divided  into  a number  of  sections  where 
the  thermal-fluid-dynamic  variables  (pressure,  velocity,  Mach 
number,  temperature  etc....)  are  computed  according  to  the  ID  flow 
hypothesis.  In  Fig.  1 a sketch  of  a reference  domain  divided  into  ID 
sections  is  presented. 


Fig.  1 : Computational  model 

Governing  equations 

The  system  of  mass,  momentum  and  energy  ID  equations,  for  each 
control  volume  (sections  i-1,  i,  i+1  in  Fig.l),  are  written  in  the 
following  conservative  form: 


d t d x 


(1) 


with 

e,  =cvT+  0.5  u2  +him, ) 

B 4mbleec!  > + Fbiade  + F friction  9 W ® ext  + Kleed  l (2) 

the  total  energy  et  is  defined  as  et  = cvT  + 0.5 w 2 . Using  the  perfect 
gas  equation  of  state: 

e,  = -3E  t + 0.5m2  P = {k~  1 )[pe,  - 0.5 pu1]  (3) 


Time-marching  technique 


In  order  to  integrate  the  nonlinear  hyperbolic  differential  system  of 
eq.l,  a time  marching  technique,  previously  developed  for  2D  and  3D 
steady/unsteady-inviscid/viscous  turbomachinery  flows  (Cravero 
1995,  Cravero- Satta  1995,  Cravero- Satta  2000),  has  been  set  up  for 
the  ID  case.  After  the  integration  in  space  the  system  of  eq.l  can  be 
written  as: 

d Qi  / \ 

-^L+(ci+A-tfi)=o  (4) 

where  the  suffix  i refers  to  the  actual  grid  node.  The  term  in 
parenthesis  is  defined  as  the  residual  Rj  and  will  be  driven  to  zero  in 
steady  flow  calculations.  The  convective  fluxes  are  integrated  with 
central  differences: 


FhlZIhl 

2Axi 


(5) 


In  order  to  prevent  odd-even  points  decoupling  and,  in  general,  to 
solve  transonic  flows  with  shock  waves,  an  artificial  dissipative  term 
is  added:  this  is  a blend  of  second  and  fourth  order  dissipative  tenns: 


A = 
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where 


5;  = B?  max(  vl+x , vt , 
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V;  = — ~ + — (9) 
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The  coefficient  02  (eq.7)  and  04  (eq.8)  are  fixed  and  their  typical 
values  are  1/4  and  1/32.  The  integration  in  time  of  eq.4  is  performed 
with  the  following  explicit  multistep  Runge-Kutta  scheme  with  N 
steps: 

ef =s(o 

fi,=ft°-alAf1*(0#) 

„ , (10) 

Of  = Qi  ) 

ad+Ao=av 

The  dissipation  term  in  the  residual  R is  computed  only  at  the  first 
step.  A four  step  scheme  with  standard  coefficients  (1/4,  1/3,  1/2,  1) 
has  been  used.  For  steady  flow  analysis  the  local  time  step  is 
computed  according  to: 

Ax 

A li  = CFL (11) 

(w  + a); 

where  the  Courant-Friedricks-Levy  has  its  limit  of  2.8  with  a four 
steps  scheme.  When  unsteady  flow's  are  computed  the  minimum  time 
step  over  the  domain  is  considered. 

Source  term  evaluation 


The  source  term  B in  eq.l  contains  the  effects  of  a mass  flowr  bleed 
(nibieed  and  hbIeed)  or  injection,  of  heat  addition/removal  (<f>ext),  of 
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channel  section  variation  (P®  /gx)>  °f  friction  at  the  walls  due  to 

viscous  effects  (Ffriction)  and  of  a turbomachinery  stage  (Fb|acle  and  W). 
Introducing  this  last  contribution  allows  the  unsteady  flow  in  a 
multistage  turbomachinery  (compressor  or  turbine)  to  be  computed 
with  many  information  of  the  flow  inside  the  machine  compared  to  a 
lumped  parameter  approach.  In  the  following  the  attention  will  be 
focused  on  the  stage  force  and  work  exchange  evaluation.  Each 
computational  station  is  identified  with  a label  that  defines  its  kind 
(i.e.  0-ducl,  1 -stator,  2-rolor).  At  present  a computational  model  for 
axial  compressor  stages  is  implemented.  When  a compressor  stage  is 
detected,  the  axial  blade  force  and  the  work  exchange  are  computed 
either  by  the  stage  characteristics  (read  from  an  external  data  file)  or 
by  a routine  implemented  into  the  code  that  returns  the  stage 
performances  in  the  actual  working  point.  The  blade  force  is  equally 
distributed  for  each  station  defining  a blade  (rotor+stator)  in  the  stage 
while  the  work  exchange  is  “spread  over”  the  stations  in  Lhe  rotor.  If 
the  internal  routine  is  called  to  compute  the  compressor  stage 
characteristics,  the  pre-stall  part  of  the  stage  characteristics  is 
computed  with  a standard  ID  approach  with  correlations  while  the 
post-stall  and  reverse  flow  parts  are  evaluated  using  the  model 
proposed  by  Moses  et  al.  (1986)  and  by  O'Brian  (1992).  A time-lag 
equation  for  the  blade  force  is  introduced  to  compute  the  actual  blade 
force  in  unsteady  calculations,  following: 


dFhl 


dt 


_ , |7  _ T7 

" ' r blade  — rblade_ steady 


(12) 


For  each  stage,  the  energy  input  to  the  overall  system  is  given  in 
terms  of  torque  coefficient.  Here  the  stage  temperature  rise 
characteristics  are  obtained  based  upon  the  overall  torque  data  and  two 
isolated  flow  points  in  rotating  stall  (all  the  stage  temperature 
characteristics  were  considered  to  be  identical). 


Fig .2:  Compression  system  geometry  - a)  Different  plenum  geometry” 
b)  section  numbering 


Inlet  and  Outlet  boundary  conditions 

At  the  inlet  section,  in  the  most  frequent  case  of  subsonic  flow,  total 
pressure  and  temperature  are  fixed  as  boundary  conditions.  The 
conservative  variables  (vector  Q of  eq.l ) are  updated  extrapolating  the 
Riemann  invariant  (Cravero,  1995)  from  the  interior  points.  At  the 
outlet  section  static  pressure  is  fixed  as  boundary  condition;  when  a 
steady  state  solution  is  required  the  conservative  variables  are 
obtained  from  the  interior  points  by  extrapolating  the  Riemann 
invariant  together  with  the  assumption  of  isentropic  flow  in  the  last 
two  sections.  When  a time  dependent  boundary  condition  is  required, 
a model  of  an  isoentropic  expansion  through  a variable  section  control 
valve  is  modeled.  The  “reduced  mass  flow”  trough  the  valve  is  used  to 
compute  the  Mach  number,  and  therefore  all  the  required  quantities: 

mvalve  = ^T=LP  = (PuA)nx  O3) 

The  total  quantities  Tt  and  pt  are  extrapolated  from  the  station  NX-1. 
The  valve  characteristic  coefficient  \x  is  computed  with  a steady  flow 
run  and  it  is,  in  general,  function  of  time  during  the  unsteady  flow 
calculation. 

APPLICATIONS 

To  test  the  model  previously  described,  the  available  geometrical 
and  experimental  data  of  the  compressor  of  Gamache  (1985)  have 
been  utilized.  The  three-stage,  low-speed  compressor  consists  of  three 
non-repeating  stages  with  a constant  cross-sectional  annulus  area  (hub 
to  tip  ratio  0.88).  From  the  cited  work  and  from  the  work  of  Eastland 
(1982),  a complete  set  of  steady  stage  pressure  characteristics  and 
corresponding  overall  steady  system  performances  are  available. 


As  already  observed  by  O'Brien  (1992),  the  referenced  three-stage 
compression  system  provides  excellent  stage  characteristic  data,  but, 
unfortunately,  overall  system  behavior  during  surge  and  rotating  stall 
was  not  available.  However,  the  rig  was  similar  to  one  used  in  a 
previous  experimental  investigation  of  surge  and  rotating  stall  by 
Greitzer  (1976).  Taking  into  account  this  aspect,  the  compression 
system  here  analyzed  has  been  configured  to  the  geometry  specified 
for  the  plenum  rig  of  Greitzer's  tests,  but  using  the  computed  stage 
characteristics  for  the  Gamache’ s compressor.  The  theoretical  stage 
performance  curves  have  been  used  for  the  following  calculations. 

For  the  analysis  of  compression  systems  Greitzer's  B parameter  was 
used  as  a reference  variable  for  the  experimental  compressor.  The 
parameter  B is  defined  as: 


(14) 


The  B value  can  be  modified  by  changing  the  plenum  volume  Yp,  as 
shown  in  Fig.2;  three  cases  were  computed:  rotating  stall  (B=0.65); 
classical  surge  (B=l  .00);  deep  surge  (B=l  .58). 

In  Fig.3  the  dynamic  responses  of  the  compression  system  with 
B=0.65-1.0-1.58  are  shown.  The  compressor  rig  throttle  was  slowly 
closed  to  the  point  of  instability  and  then  held  constant.  In  case  of 
B=0.65  the  system  became  unstable  at  the  uniform  flow  stall  point, 
and  then  traversed  to  rotating  stall.  A second  calculation  has  been 
done  with  the  compressor  rig  reconfigured  to  operate  at  B=1.00  (see 
Fig,2a).  In  this  condition  the  system  exhibited  surge  cycles.  In  Fig.3 
the  time-dependent  distributions  of  static  pressure  in  different  stations 
along  the  machine  clearly  show  the  different  wave  amplitude  and  the 
need  for  a distributed  model  to  simulate  the  system.  Another  test  has 
been  considered  with  the  compressor  rig  reconfigured  to  operate  at 
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B=1.58  (see  Fig.2a).  In  this  configuration  the  system  exhibited  deep 
surge  cycles. 

B=0.65 


B=1.58 


Fig. 3:  Pressure  and  velocity  variations  in  different  sections  in  case  of 
B=0.65  - B=1 .0  - B=1.58 

The  large  amplitude  variations  of  static  pressure  inside  the  compressor 
can  be  detected  by  the  ID  code  as  shown  in  Fig. 3.  It  is  clear  from  the 
examples  that  the  present  ID  model,  compared  to  a lumped  model 
approach,  can  give  detailed  information  about  the  flow  behavior  inside 
the  compressor  (i.e.  interstage  fluid  dynamic  data)  and  inside  the 
plenum.  This  is  an  important  feature  of  the  model  that  allows  the 
influence  of  different  active  control  strategies  (interstage  bleed, 
cooling  mass  flow  rate,  valves  ....)  to  be  analyzed  and  discussed.  In 
Fig.4  the  time  dependent  velocity  distributions  within  the  system, 
obtained  with  different  configurations  (B-0. 65- 1.0- 1.58),  are  shown 
in  a 3D  plot  to  highlight  the  strong  variations  inside  the  compression 
system.  This  kind  of  inspections  can  be  used  to  understand  the 
optimum  placing  of  transducers  and  active  control  devices. 


Fig.4:  3D  views  of  velocity  distributions  with  respect  to  time  and 
space-  a)  B=0.65  - b)  B^l.O  - c)  B=1.58 
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ACTIVE  CONTROL  STRATEGY 

Compression  systems  have  their  operating  range  limited  at  high  flow 
rate  by  choking  and  at  low  flow  rates  by  occurrence  of  aerodynamic 
instabilities,  namely  rotating  stall  and/or  surge.  If  we  concentrate  our 
attention  on  a cascade:  as  flow  rate  decreases,  the  incidence  angle 
grows  up  until  flow  separates  so  that  the  cascade  pressure  rise  and 
efficiency  fall  out.  Depending  on  the  compressor  and  plant 
characteristics  there  will  be  stall  or  surge.  The  former  is  a non- 
axisymmetric  perturbations  with  one  or  more  cells  with  lower 
meridional  velocity  and  “clean”  flow  zones,  but  stationary 
circumferentially  averaged  flow  rate;  virtually  it  might  be  not 
dangerous.  On  the  contrary,  the  latter  involves  strong  perturbations  in 
flow  rate,  sometimes  exhibiting  reverse  flow  periods,  and  it  can 
seriously  damage  both  compressor  and  plant.  Compressor  map  and 
expander  line  (if  it  is  a throttle  valve)  or  map  (if  it  is  a turbine),  but 
also  volumes  (tubes,  burner  and  plenum)  between  them  play  a role  to 
lead  to  stall  or  surge.  Greitzer  (1976)  pointed  out  that  it  is  not  only  a 
“local”  matter,  so  that  surge  instability  can  arise  even  if  the  slopes  of 
the  characteristics  seem  “good”.  Despite  of  this,  most  researchers 
utilize  linear  models  (Ffowcs  Williams  and  Huang,  1989;  Gysling  et 
al.,  1991;  Pinsley  et  al.,  1991;  Simon  et  al.,  1993;  Jungowsky  et  al., 
1994).  However,  to  avoid  surge,  traditional  devices  keep  the 
compression  systems  away  to  work  close  to  low  incidence  angles,  i.e., 
low  flow  rates.  By  the  way,  such  zones  would  exhibit  very  good 
pressure  rise  and  efficiency. 

In  the  last  decade  new  techniques  have  been  searched,  which  are 
more  effective  than  purely  inhibiting  dangerous  zones.  Formerly 
Greitzer  found  that  instability  was  related  to  a rate  of  energy,  so  that 
many  researchers  thought  that,  if  this  energy  could  be  dissipated, 
instability  had  to  disappear. 

Two  families  of  techniques  were  suggested:  passive  and  active 
control.  The  former  is  based  on  the  aeroelastic  coupling  of  the 
compression  system  with  a mass-spring-damper  device  (Gysling  et  al., 
1991)  or  a hydraulic  oscillator  (Amulfi  et  al.,  2000),  able  to  catch  and 
dissipate  the  instability  energy.  Unfortunately,  this  approach  involves 
serious  practical  problems  and  actually  no  industrial  plant  has  been 
equipped  with  such  a device.  The  latter  family  is  based  on  a classic 
feedback  control  system,  i.e.,  a sensor/'actuator  pair  and  a suitable 
control  law.  Ffowcs  Williams  and  Huang  (1989)  stabilized  the 
unstable  flow  in  a small  compression  system  by  using  a pressure 
sensor  located  in  the  plenum  and  a movable  wall  as  an  actuator  (a  sort 
of  loudspeaker  diaphragm).  One  hardly  can  think  of  an  industrial  size 
compression  unit  with  a device  like  this.  Pinsley  et  al.  (1991) 
suggested  a more  promising  control  system,  where  the  plenum 
pressure  signal  modulates  the  compressor  exit  area  by  a throttle  valve, 
located  downstream  of  the  plenum.  Although  they  experimentally 
validated  their  idea  with  a small  centrifugal  compressor  for 
automotive  turbo-charging,  it  seems  to  be  suitable  to  large  industrial 
units  too.  Many  other  researchers  concentrate  their  study  on  active 
stabilization  of  surge  in  small  units  (for  example,  Jungowsky  et  al., 
1994,  Komatsubara  and  Mizuki,  1995). 

As  far  as  axial  compressors,  and  particularly  aero-engine  units,  there 
is  a lack  of  experimental  work.  Cargill  and  Freeman  (1991)  describe 
an  active  control  scheme.  The  control  effectiveness  needs  two 
requirements:  very  fast  response  to  frequency  and  a large  number  of 
sensors  and  actuators.  The  former  conflicts  with  the  use  of  analogic  or 
digital  filters:  they  have  to  be  utilized  in  order  to  remove  both  noise 
and  DC  component  of  the  experimental  signals,  but  cause  some  phase 
shift.  The  latter  makes  the  plant  complex  and  hardly  manageable. 
However  Montazeri  et  al.  (1996)  wrote  a theoretical  paper  on  an 


active  control  device  using  one  actuator  only.  It  deals  with  an  air 
bleeding  strategy  and  the  location  of  the  bleeder  was  optimized. 

Several  control  strategies  can  be  considered  for  the  active 
stabilization  of  a compressor  system  under  transient  operations  that 
can  cause  the  compressor  instability.  With  the  computational  model 
described  above  the  effect  of  different  interventions  (distributed  mass 
flow  bleed,  variable  geometry  devices,  control  valves  etc. . .)  made  to 
stabilize  the  compressor  system  can  be  evaluated  and  the  most 
efficient  active  control  strategy  can  be  outlined. 
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Fig.5:  Compression  system  with  control  valve  and  main  valve  with  its 
operational  law 

In  this  paper  a preliminary  application  of  a simple  active  control 
strategy  has  been  considered  in  order  to  verify  the  stability  properties 
of  the  iterative  solution  method  and  to  estimate  the  potential 
application  of  the  code  to  more  complex  control  devices.  The 
compression  system  operation  point  is  set  by  the  main  valve  (device 
D1  of  Fig.5)  whose  operational  law  is  defined  by  the  initial  (Al)  and 
final  (A2)  values  of  the  channel  section  (related  to  the  coefficient  p of 
eq.13)  and  by  the  time  step  dt  of  the  operation. 
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Fig.  6:  Compression  system  instability,  a)  flow  coefficient,  b) 
static  pressure 

A control  valve  (device  D2  in  Fig.5)  is  introduced  in  the  exit  duct  in 
order  to  stabilize,  with  a defined  control  law,  the  system.  The 
following  simple  control  law  for  the  valve  B has  been  introduced: 

AA  = Kdp/dt  (14) 

where  AA  is  the  section  variation  and  the  pressure  gradient  is  taken  at 
the  compressor  exit. 
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Fig. 7:  Compressor  system  stabilization,  a)  flow  coefficient,  b)  static 
pressure,  c)  control  valve  operation 


If,  starting  from  the  nominal  operating  condition,  a reduction  of 
about  25%  in  the  exit  duct  section  is  operated  in  0.1  s by  the  main 
valve  Dl,  the  compression  system  goes  into  a surge  behavior.  This  is 
shown  in  Fig.6  where  the  flow  coefficient  and  the  static  pressure 
(referred  to  the  starting  values)  transient  behavior  are  plotted  for  the 
inlet  (solid  line)  and  outlet  (dashed  line)  compressor  stations.  With  the 
same  transient  operation  set  by  the  main  valve  D l the  compression 
system  is  stabilized  by  acting  on  the  control  valve  D2  with  the  law  set 
by  eq.14.  In  Fig.7a-7b  the  flow  coefficient  and  the  static  pressure  are 
plotted  against  the  time  for  the  inlet  (solid  line)  and  outlet  (dashed 
line)  compressor  stations.  The  initial  instability  is  damped  after  a short 
time  by  the  active  intervention  of  the  control  valve  D2.  The  control 
valve  is  activated  after  the  transient  operation  as  shown  in  Fig. 7c 
where  the  unsteady  section  variation  is  plotted.  As  predicted  by  the 
model  very  small  variations  are  required,  but  with  very  high  gradients. 

This  is  a preliminary  simulation  for  an  active  control  strategy  using 
the  distributed  computational  model  and  the  calculation  has  been 
performed  mainly  to  verify  the  stability  of  the  explicit  scheme  when 
introducing  variable  sections.  Actually  no  attention  has  been  devoted 
to  the  practical  aspects  for  the  implementation  of  the  above  control 
strategy  in  a real  system.  In  future  work  the  effect  of  different  control 
strategies  on  the  active  stabilization  of  the  compression  system  will  be 
assessed  with  respect  to  their  practical  implementation. 


CONCLUSIONS 


A time-marching  one  dimensional  simulation  model  has  been  set  up 
for  the  unsteady  aerodynamic  analysis  of  a gas  turbine  compressor. 
The  model  considers  the  system  with  its  real  geometrical 
configuration  and  it  is  therefore  suitable  to  the  computational  analysis 
of  the  effect  of  different  active  control  devices  on  the  compressor 
performances.  After  the  verification  of  the  ability  of  the  code  to  give 
reliable  solutions  for  a multistage  axial  compressor,  even  with  an 
unstable  operating  mode,  a preliminary  active  control  strategy  has 
been  implemented  in  order  to  check  for  the  stability  of  the  numerical 


scheme.  The  code  has  demonstrated  good  stability  properties  and  it 
will  be  therefore  used  in  future  works  to  understand  the  effectiveness 
of  different  active  control  strategies.  The  fluid-dynamic  solver  has 
been  written  as  a general  solution  tool  for  one  dimensional  problems 
in  order  to  be  able  to  extend  the  analysis  to  a complete  gas  turbine 
engine  with  minor  modifications  to  the  code.  This  approach  will  be 
used  to  simulate  the  effect  of  different  fluid-dynamic  instabilities 
(compressor  surge,  burner  humming  etc...)  on  the  whole  system  with 
allowance  of  active  control  strategies, 
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